Dissipative Taylor-Couette flows under the influence of helical magnetic fields 
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The linear stability of MHD Taylor-Couette flows in axially unbounded cylinders is considered, for magnetic 
Prandtl number unity. Magnetic fields varying from purely axial to purely azimuthal are imposed, with a general 
helical field parameterized by /3 = / B z . We map out the transition from the standard MRI for j3 = to 
the nonaxisymmetric Azimuthal MagnetoRotational Instability (AMRI) for f3 — > oo. For finite /?, positive and 
negative wave numbers m, corresponding to right and left spirals, are no longer identical. The transition from 
/3 — to ft — > oo includes all the possible forms of MRI with axisymmetric and nonaxisymmetric modes. For 
the nonaxisymmetric modes, the most unstable mode spirals in the opposite direction to the background field. 
The standard (/? = 0) MRI is axisymmetric for weak fields (including the instability with the lowest Reynolds 
number) but is nonaxisymmetric for stronger fields. 

If the azimuthal field is due in part to an axial current flowing through the fluid itself (and not just along the 
central axis), then it is also unstable to the nonaxisymmetric Tayler instability, which is most effective without 
rotation. For large /3 this instability has wavenumber m = 1, whereas for /3 ~ 1 m = 2 is most unstable. The 
most unstable mode spirals in the same direction as the background field. 

PACS numbers: 47.20.Ft, 47.20.-k, 47.65.+a 



I. INTRODUCTION 

The longstanding problem of the generation of turbulence 
in various hydrodynamically stable situations has found a so- 
lution in recent years with the MHD shear flow instability, the 
so-called magnetorotational instability (MRI), in which the 
presence of a magnetic field has a destabilizing effect on a 
differentially rotating flow with outward decreasing angular 
velocity but increasing angular momentum |l|,|2|]. 

In the absence of MHD effects, according to the Rayleigh 
criterion, an ideal flow is stable against axisymmetric pertur- 
bations whenever the specific angular momentum increases 
outward 

^(i? 2 tt) 2 > 0, (1) 

where SI is the angular velocity, and (R, <fr, z) are cylindrical 
coordinates. In the presence of an azimuthal magnetic field 
Ba,, this criterion is modified to 
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where p,o is the permeability and p the density . Note also 
that this criterion is both necessary and sufficient for (axisym- 
metric) stability. In particular, all ideal flows can thus be 
destabilized, by azimuthal magnetic fields with the right pro- 
files and amplitudes. 

On the other hand, for nonaxisymmetric modes, one has 
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as the necessary and sufficient condition for stability of an 
ideal fluid at rest @], Outwardly increasing fields are there- 
fore unstable. If (01 is violated, the most unstable mode has 
azimuthal wavenumber m = 1 . 

The rich variety of nonaxisymmetric instabilities can be 
demonstrated by the addition of a differential rotation. In this 
case even the current-free (within the fluid) profile ocl/R 
(which according to (O is stable for 57 = 0) can become un- 
stable. Even for a differential rotation that by itself would be 
stable according to [T] the combination of £7 and B^ oc 1/R 
can be unstable to to = 1 perturbations (see Fig. [2j. We have 
called this phenomenon the Azimuthal MagnetoRotational In- 
stability (AMRI). It has even been demonstrated that it should 
be possible to observe the AMRI in laboratory experiments, 

B. 

Further new phenomena appear if an axial field is added, 
yielding a spiral, or helical total field. In this case only a suf- 
ficient condition for stability against axisymmetric perturba- 
tions is known. In the absence of rotation this is 



(4) 



(see also Eq. ©). Including rotation, this was extended |0] to 
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For the current-free field B^ oc 1 /R, only superrotating flows 
with dfl/dR > are stable. Indeed, we have demonstrated 
that dissipative Taylor-Couette flows beyond the Rayleigh 
limit for centrifugal instability can easily be destabilized by 
helical magnetic fields with such a current-free azimuthal 
component [7]. The resulting axisymmetric traveling wave in- 
stability has become known as the Helical MagnetoRotational 
Instability (HMRI), and has been obtained in the PROMISE 
experiment J8l|9|]. 
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In the PROMISE experiment the azimuthal field is oc 
1 /R. In this paper we will also consider the generalization to 
Bff, — asR + bs I R, where the extra term asR corresponds 
to an axial electric current running through the fluid as well, 
and hence opens the possibility of current-induced (Tayler) in- 
stabilities. The resulting (nonaxisymmetric) instabilities may 
also be modified by adding either a differential rotation or an 
axial magnetic field. 

One might suppose that adding an axial field would be im- 
portant only if its amplitude is of the same order as that of the 
azimuthal field. Chandrasekhar ifloll showed that for il = 0, 
a sufficiently strong axial field will always suppress any ax- 
isymmetric instabilities of an azimuthal field, by deriving the 
stability condition 



IBl > 



&±(&Bl) dR, 



(6) 



where I > and £r is the (purely real) radial eigenfunc- 
tions. (Note how (O reduces to <(4j for B z = 0.) However, we 
will show that the influence of B z cannot be ignored even for 
rather small values. 

We will find that, depending on the magnitudes of the im- 
posed differential rotation and magnetic fields, the field may 
either stabilize or destabilize the differential rotation, and the 
most unstable mode may be either the axisymmetric Taylor 
vortex flow (the SMRI or HMRI), or the nonaxisymmetric 
AMRI, or the nonaxisymmetric Tayler instability. In com- 
bined axial and azimuthal fields, we will also show that the 
nonaxisymmetric modes differ between m and — m, corre- 
sponding to left and right spirals. As first pointed out by Hill , 
if the imposed field has both axial and azimuthal components, 
the system no longer exhibits ±z symmetry. For axisymmet- 
ric modes, the consequence of this is that what were previ- 
ously stationary modes (SMRI) become oscillatory, traveling 
wave modes (HMRI). For nonaxisymmetric modes, breaking 
the ±z symmetry of the basic state breaks the ±m symmetry 
of the instabilities. Physically this corresponds to the fact that 
modes spiraling either in the same or the opposite sense to the 
spiral structure of the basic state are indeed different. This 
±m symmetry breaking is also a convenient distinguishing 
feature between the AMRI and the Tayler instabilities; for the 
AMRI the most unstable mode spirals in the opposite sense 
to the imposed field, for the Tayler instabilities in the same 
sense. 

Finally, in order to produce benchmarks for the application 
of incompressible 3D MHD codes, in this work we will focus 
primarily on magnetic Prandtl number Pm = 1. 



Developing the disturbances into normal modes, the solutions 
of the linearized MHD equations are considered in the form 



F = F(R)exp(i(kz + mcf) + uit)), 



(8) 



where F is any of the velocity, pressure, or magnetic field 
disturbances. 

The governing equations are 

OIL \. 

— + {U -V)u + (u-V)U = --Vp + vAu+ (9) 

H — —curl bx B + — !— curl B x b, (10) 
PoP PoP 



db 

dt 



and 



curl(u x B) + curl(J7 x b) + rjAb, (11) 



div it = div b = 0, (12) 



where u is the perturbed velocity, b the magnetic field, p the 
pressure perturbation, v is the kinematic viscosity and r\ the 
magnetic diffusivity. 
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FIG. 1 : The basic state azimuthal field with prescribed values at the 
inner and outer cylinders. 



The stationary background solution is 

fi = an + -j^, B^ = asR + 
where an, bn, qb and bs are constants defined by 



(13) 



II. THE EQUATIONS 

We are interested in the linear stability of the background 
field B = (0,B^(R),B ), with B =const, and the flow 
U = (0, Rft(R), 0). The perturbed state of the system is 
described by 



ur, M0, u z , p, b R , &0, b z 



(7) 
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i?i n and i?out are the radii of the inner and outer cylinders, 
f2; n and f2 ou t are their rotation rates, and B ln and B out the az- 
imuthal magnetic fields at the inner and outer cylinders. The 
possible magnetic field solutions are plotted in Fig. [T] Note 
that - unlike f2, where Q- m and £l ou t are the physically rele- 
vant quantities - for the fundamental quantities are not so 
much i?i n and B out , but rather a B and b B themselves. In par- 
ticular, a field of the form b B /R is generated by running an 
axial current only through the inner region R < R- m , whereas 
a field of the form clbR is generated by running a uniform ax- 
ial current through the entire region R < i? ou t, including the 
fluid. 

Given the z-component of the electric current, curl z £? = 
2clb, one finds for the current helicity of the background field 



curLB • B = 2a B B , 



(16) 



which may be either positive or negative (and of course van- 
ishes for the current-free case a B = 0). However, both signs 
yield the same instability curves, merely with the previously 
mentioned left and right spirals interchanged. 

The inner value B ln is normalized with the uniform vertical 
field, i.e. 



/3 = 



-Bin 
Bn 



For i^b = 1 we have 



curl-B • B = 



2£Bl 

3 Ri n 



(17) 



(18) 



The sign of f3 thus determines the sign of the helicity of the 
background field. And again, interchanging ±/3 simply inter- 
changes left and right spirals ±m. 

As usual, the toroidal field amplitude is measured by the 
Hartmann number 



Ha = 



B iD R 



(19) 



Rq = J Bin (-Rout ^ -Rin) is used as the unit of length, tj/Rq 
as the unit of velocity and B- m as the unit of the azimuthal 
fields. Frequencies, including the rotation f2, are normalized 
with the inner rotation rate f2i n . The magnetic Reynolds num- 
ber Rm is defined as 



Rm = 



$2; n Rq 



(20) 



The Lundquist number S is defined by S = Ha ■ \/Pm. 

The boundary conditions associated with the perturbation 
equations are no-slip for u, 



UR = M0 = u z = 0, 
and perfectly conducting for b, 

db^/dR + b^/R = b R = 0. 



(21) 



(22) 



These boundary conditions hold for both R = R{ n and R = 

Rout- 



If we consider a constant phase of a nonaxisymmetric pat- 
tern, Eq. dHJ yields 
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(23) 



The first relation describes the phase velocity of the modes in 
the axial direction, the second in the azimuthal direction (and 
only exists for nonaxisymmetric modes). Obviously the wave 
is traveling upwards if the real part of the eigenfrequency is 
negative. 

At a fixed time the phase relations (l23l can also be written 

as 



dz/d(j) = —m/k. 



(24) 



Now, k and m are both real numbers, and without loss of gen- 
erality one of them can be taken to be positive, say k. The 
other one, m in this case, must be allowed to have both signs 
though. Negative m describe right-hand spirals (marked here 
by R), and positive m describe left-hand spirals (marked here 
by L). If the axisymmetric background field possesses positive 
B z and B^ (as used for the calculations here) then its current 
helicity is positive, or equivalently, it forms a right spiral. 



HI. FROM AMRI TO HMRI 

We begin with a purely azimuthal field, and no electric cur- 
rents within the fluid, that is, B$ <xl/R. Figure|2]presents re- 
sults for (Uq = 0.5, showing that for Ha > 100 and Re > 200 
there exists an m = 1 nonaxisymmetric instability. Note also 
how both the upper and lower branches of the instability curve 
tilt to the right, that is, have a positive slope dRe / dHa. For 
a given Hartmann number, the instability therefore only ex- 
ists within a finite range of Reynolds numbers. If Re is too 
large, the instability disappears again as a consequence of the 
suppressing action that differential rotation often has on non- 
axisymmetric modes. 

We consider next a purely axial field, the so-called stan- 
dard MRI. In this case both axisymmetric and nonaxisymmet- 
ric modes may be excited, with the axisymmetric mode being 
the one with the overall lowest Reynolds number (Fig. |3). 
For Pm = 1 this overall minimum occurs for Ha ~ 10 and 
Re ~ 80. However, for sufficiently large Ha the m = 1 non- 
axisymmetric mode is actually preferred over the axisymmet- 
ric mode. The standard MRI for purely axial fields is there- 
fore not necessarily an axisymmetric mode. The axisymmet- 
ric mode only dominates for sufficiently weak fields, includ- 
ing also the global minimum Re value. It also dominates the 
entire weak-field branch of the instability curves (Pig. [3j. The 
axisymmetric mode here tilts to the left, whereas the nonax- 
isymmetric mode tilts to the right, as before in Fig. [2] 

Figure [4] finally shows results combining azimuthal and ax- 
ial fields, focusing in this case on (3 = 2 (so a right-handed 
helicity). We see the same general pattern as before: only the 
weak-field branch of the m = mode tilts to the left; all non- 
axisymmetric modes tilt to the right. Up to Ha rj 50 the ax- 
isymmetric mode is preferred, just as before for the standard 
MRI. For Ha > 50 the m = 1 right spiral is preferred. 
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FIG. 2: Azimuthal MagnetoRotational Instability (AMRI): The in- 
stability curve for current-free azimuthal magnetic fields. The most 
unstable mode is nonaxisymmetric, with m — 1. Axisymmetric 
instabilities do not exist at all in this configuration. \ib = 0.5, 
fin = 0.5. Pm = 1. 
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FIG. 3: Standard MRI: The instability curves for uniform axial mag- 
netic fields. For Ha < 75 the axisymmetric mode is preferred; for 
Ha > 75 the m = 1 nonaxisymmetric mode is preferred. The 
global minimum Re value is for the axisymmetric mode, = 0.5, 
Pm = 1. 



Figure|4]also demonstrates that the (axisymmetric) standard 
MRI and the (nonaxisymmetric) AMRI are the basic elements 
which both appear, with different weights, if the background 
field has a spiral geometry. From this point of view instabil- 
ities in helical fields are simply a mixture of these two ba- 
sic elements. More specifically, one finds that the weak-field 
branch of the instability in Fig. [4] is very similar to the weak- 
field branch of the standard MRI (Fig. [3) while the strong-field 
branch strongly resembles the strong-field branch of AMRI 
(Fig. |2j. The minimum is always obtained for the axisym- 
metric mode of the standard MRI. The only difference be- 
tween the standard MRI and HMRI is the different character 
of the eigenfrequencies: the SMRI is stationary, whereas the 
HMRI is oscillatory, as a necessary consequence of the ±z 



symmetry -breaking fiUl . It is precisely this oscillatory nature 
of the HMRI that has been used to identify it in the PROMISE 
experiment (fa |9|] . 
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FIG. 4: Helical magnetorotational instability (HMRI): The instabil- 
ity curves for current-free (/i_g = 0.5) helical fields with ft = 2. 
(j,n = 0.5. The dotted line gives the axisymmetric mode. The solid 
lines are marked with the mode number m and the type of helicity (L 
left, R right). Note that for stronger fields the most unstable mode is 
nonaxisymmetric with the same helicity as the helicity of the back- 
ground field. The mode with the lowest Re is always axisymmetric. 
Top: Pm = 1, bottom: Pm = 0.01. 

Note finally that taking Pm = 1 greatly simplifies the re- 
sults, and indeed eliminates some particularly interesting re- 
sults. As we have previously demonstrated, both the (axisym- 
metric) HMRI 10] as well as the (nonaxisymmetric) AMRI |5| 
have the property that their scalings with Pm vary dramatically 
with fiQ. For /j,q only somewhat greater than the Rayleigh 
value, both modes have Ha and Re as the relevant measures 
of field strength an d rot ation rates, whereas for greater values 
of /in, S = Ha • y/Pm and Rm = Re • Pm are the relevant 
measures. For small Pm the differences can thus be huge. 
Insulating versus conducting boundaries can also have a sur- 
prisingly large influence on this transition from one scaling to 
another d. 
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IV. FIELDS WITH CURRENT-HELICITY 



A. Steep rotation law 
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FIG. 5: Tayler instability (TI) of a toroidal field under the influence 
of differential rotation with stationary outer cylinder. The solid curve 
is m — 1, the dotted curve is m — 0. fiB = 1, fin = 0. Pm = 1. 

We begin by considering the stability of purely toroidal 
fields, and differential rotation profiles with a stationary outer 
cylinder. There are then three classical results known: First, in 
the absence of any fields, axisymmetric Taylor vortices arise 
at Re = 68, and nonaxisymmetric instabilities at Re = 75. 
Second, in the absence of any rotation, m = 1 Tayler in- 
stabilities arise at Ha = 150. Figure [5] shows how these 
results are linked when both Ha and Re are non-zero. For 
Ha very small, the axisymmetric Taylor vortex mode is stabi- 
lized, whereas the nonaxisymmetric mode is eventually desta- 
bilized, and connects smoothly to the pure Tayler instability. 

The next step is to add a uniform axial field to the azimuthal 
field, with (say) positive polarity. The background field then 
has a positive helicity, that is, it spirals to the right. If the axial 
field is weak, e.g. with j3 — 100, then the marginal instability 
curves (Fig. |6l top) strongly resemble the map for f3 — > oo 
(Fig. [5]). The main differences are i) the slightly smaller Hart- 
mann number of the toroidal field, and ii) the splitting of the 
spiral modes m = 1 and m = — 1 into two curves with differ- 
ent helicity (R and L). The left-hand modes require a greater 
rotation than the right-hand modes. For background fields 
with positive helicity, we thus find that the right spirals are 
preferred, whereas for background fields with negative helic- 
ity, R and L would be exchanged, and the left spirals would 
be preferred. 

For j3 — 10 the differences between the L and R modes 
for given m increase, so that the pure Tayler instability exists 
only as the 1R mode. The 1L mode no longer connects to 
Re = 0, and is not the most unstable mode anywhere in the 
given domain (Fig. |6l middle). 

The 1R mode also dominates for j3 of order unity. There 
is, however, an interesting particularity in this case. For very 
slow rotation, a 2R mode reduces the stability domain. For 
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FIG. 6: Instability curves for magnetic fields with positive current 
helicity. The curves are marked with their azimuthal wave num- 
bers m. The shaded areas are the stable regions. Without magnetic 
fields the curves always start at Re = 68 for the m = Taylor vor- 
tices. Note the dominance of the axisymmetric modes (dotted) also 
for weak magnetic fields, unless for very high /3 the field becomes 
nearly toroidal. Pm = 1, fiB = 1, fin = 0. 



Re ~ 0, and in a limited range of Ha (Ha ~ 100... 130), this 
mode forms the first instability (see [13]). A small amount of 
differential rotation, however, brings the system back to the 
1R instability. 

For models with helical fields and steep rotation laws (with 
stationary outer cylinder), we indeed find the expected split- 
ting between right and left spiral instabilities. If the axisym- 
metric background field is right-handed, then the first unstable 
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mode is also right-handed. The corresponding critical mag- 
netic field strength is reduced compared to the TI of purely 
toroidal fields. If both magnetic field components are of the 
same order then the 2R mode is found to destabilize the sys- 
tem at the strong-field side of the stability domain, but only 
for very slow rotation. The differential rotation basically lim- 
its the action of this 2R mode. 



B. Flat rotation law 

For weak magnetic fields and the steep rotation law, the ax- 
isymmetric Taylor vortex mode is the most easily excited in- 
stability. For a sufficiently flat rotation law the non-magnetic 
Taylor vortices necessarily disappear, and a critical Reynolds 
number no longer exists for Ha = 0. For the flat rotation 
law with fiQ = 0.5, and the nearly uniform toroidal field with 
Hb = 1, the instability curves for purely toroidal fields are 
given by Fig. [7] Both of the previous instabilities appear in 
this case: TI exists in the lower right corner, and AMRI ex- 
ists in the upper left corner. The AMRI arises from the term 
bs/R in the magnetic background field profile ( TT3T >. while the 
current-driven TI is due to the term a^i?. The two instabili- 
ties are separated by a stable branch with Re w Ha, where the 
differential rotation stabilizes the TI. 
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FIG. 7: Flat rotation law (/j,q = 0.5) and B z = 0. The curves are 
marked with their mode number m. The AMRI occurs in the upper 
left comer, and the TI in the lower right corner, /is = 1, Pm = 1 

For /3 = oo the modes with positive and negative m are 
degenerate. At the weak-field limit the line for m = 2 even 
crosses the line for m = 1. Nevertheless, the AMRI solu- 
tion with the lowest Reynolds number is a nonaxisymmetric 
mode with m = 1. We find that this remains true for heli- 
cal background fields with large j3, but for j3 of order unity 
and smaller the m = mode yields the instability with the 
lowest Reynolds number (Fig. [8j — as is also true for the 
standard MRI and HMRI. The transition from nonaxisymme- 
try to axisymmetry can be accomplished simply by increasing 
the axial component of the background field. It is thus clear 
that there is a smooth transition from one form of the MRI in 



TC flows to the next. The same is true for the corresponding 
eigenfrequencies, which develop from real values (for stan- 
dard MRI) to complex values (in all other cases). 

On the other hand, if large-scale electric currents flow 
through the fluid, a critical Hartmann number exists for Re = 
0, similar to Fig. [5] where the system is also unstable even 
for Re = 0. In this case the critical Hartmann number is un- 
changed; it is again Ha = 150 for purely toroidal fields, i.e. 
[3 = oo (Figs. |5J |7J. This value does not even depend on 
the magnetic Prandtl number. For increasing j3, however, the 
critical Hartmann number is reduced to about 100. The most 
unstable mode is 1R for f3 > 10, but is again 2R for j3 of 
order unity. This result holds for very weak differential rota- 
tion; only then a mode higher than m = 1 plays a role in the 
transition from stability to instability. 

For background fields with positive helicity the TI favors 
instability patterns with right spirals. 

The instability curves of the weak-field, or diffusion- 
dominated (AMRI) limit also show a characteristic behavior. 
For large j3 it is formed by the nonaxisymmetric modes, while 
for small j3 the axisymmetric mode prevails. Consequently, 
the slopes of the lines change from positive for the nonaxisym- 
metric modes to negative for the axisymmetric modes (Fig. [8]). 
Again, the transition from AMRI to standard MRI becomes 
clear by variation of f3. If the preferred modes are nonaxisym- 
metric (for large p), then the spirals are always left-handed. 
The different mode pattern is the characteristic difference to 
the preferred modes in the TI domain. 



C. No rotation 

In general, for given Hartmann number the differential rota- 
tion stabilizes the Tayler instability which also exists without 
any rotation. On the other hand, we have shown that the crit- 
ical Hartmann numbers for nonrotating containers do not de- 
pend on the given value of the magnetic Prandtl number Pm, 
11411 . Hence, the results given in Fig. [9] for Re = and for 
m = — 1 ... — 5 are also valid for the small magnetic Prandtl 
numbers of liquid metals such as sodium or gallium which are 
used in the laboratory. 

The question about the critical Hartmann numbers for /3 < 1 
arises if the azimuthal mode number m is varied. Generally 
the mode with m = —2 dominates but for j3 < 0.4 the mode 
with m — — 3 starts to be preferred. It may happen that even 
higher m appear to be preferred for even smaller j3. However, 
it will only happen for so high values of the Hartmann number 
( > 2500) that i) laboratory experiments are impossible and 
ii) numerical investigations with differential rotation included 
- which in particular stabilizes higher m - are not possible. 
The basic result of the calculations is that the reduction of 
the increase of the axial field component (ft < 1) acts strongly 
stabilizing. This the more as the normalized differences of 
the critical Hartmann numbers for various m become smaller 
and smaller. These results do not change if formulated with 
the Hartmann number of the axial field rather than with the 
Hartmann number of the toroidal field. The total energy which 
is necessary to excite TI strongly grows with decreasing j3. 
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FIG. 8: The same as in Fig.|7]but for finite values of fi: Left: /3 = 100 (top) and /3 = 10 (bottom), Right: /3 = 2 (top) and /3 = 0.1 (bottom). 
The curves are marked with the azimuthal mode numbers m, the curves for m — are dotted. The notation R (right spiral) stands for negative 
m and the notation L (left spiral) stands for positive m, It is /i_g = 1, Pm = 1, 



Absolutely no instability remains for the limit /3 — > 0. 

We know from previous calculations that for £1 = an al- 
most homogenous toroidal field (/is = 1) becomes unsta- 
ble against disturbances with azimuthal number m = — 1 for 
Ha > 150. If an axial field is added then the critical Hartmann 
number is reduced, i.e. the toroidal field is destabilized by the 
axial component. While for B z = no preferred helicity 
exists for the instability pattern with axial field the resulting 
spiral geometry is the same as that of the background field. 
We also find that a total minimum of the critical Hartmann 
number exists for (3 < 10 (typical values of the experiment 
PROMISE) where the mode with m = 1 is the most unstable 
one. If the axial field starts to dominate for f3 < 2 then the 
critical Hartmann numbers are growing, i.e. system becomes 
more and more stable (see Fig. |9j. 

V. NONLINEAR SIMULATIONS 

The previous results have all been purely linear onset cal- 
culations, in which the governing equations are reduced to a 
linear, one-dimensional eigenvalue problem. It is also of in- 
terest to study the nonlinear equilibration of some of these 



modes, which we do with a three-dimensional spectral MHD 
code IU5I1 . The code is based on Fourier modes in 0; for each 
Fourier mode the (R, z) structure is discretized by standard 
spectral element methods involving Legendre polynomials 
01611 . For Reynolds numbers only slightly beyond the linear 
onset, the solutions do not develop much structure yet, so 16 
Fourier modes were sufficient in <f>. For the (R, z) structure, 
we typically used 2 spectral elements in R, and 14 in z (where 
periodicity with a domain height of V = 27r(i? out — R[ n ) was 
enforced), with a polynomial order between 12 and 18. The 
time-stepping uses third-order Adams-Bashforth for the non- 
linear terms, and second-order Crank-Nicolson for the diffu- 
sive terms. Boundary conditions are as before, no-slip for U, 
and perfectly conducting for B. Initial conditions are the ba- 
sic Couette profile for U, and random perturbations of size 
10~ 6 .B in for B. 

We begin by verifying that transforming (3 — >• — f3 has the 
expected result. Positive/negative (3 do indeed yield right/left 
mirror image spirals, verifying the linear onset conclusion that 
these instabilities spiral in the same sense as the imposed field. 

The simulations concern the linear onset curves for flat ro- 
tation law (see Fig. [8) both for AMRI and TI. Two examples 
are given for each instability, to probe the spatial pattern and 
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FIG. 9: The critical Hartmann number | |19I > for m = —1 ... — 5 
and /3 £ 10. An increasing dominance of the axial magnetic field 
component acts stabilizing. The curves do not depend on Pm. /is = 
1. 

the resulting field strength of the modes. The helical structure 
of all solutions is clearly visible, dominated by low Fourier 
modes m — 1 and/or m = 2 in agreement with the linear 
analysis. The solutions are stationary, except for a drift in the 
azimuthal direction. 

Figure [TO] concerns the AMRI domain. The value of (3 — 



10 is fixed, but the location in the instability diagram differs 
slightly. The top row shows the AMRI just for the minimum in 
Fig. [8] while the bottom row shows higher parameter values. 
In both cases though we see the expected m = 1 left spirals, 
in agreement with the linear results. 

The simulations lead to a further basic result. By consid- 
ering the maximum values of the radial and azimuthal com- 
ponents, a distinct anticorrelation becomes visible. The az- 
imuthal component has its maximum where the radial com- 
ponent has its minimum. The azimuthal average of bsb^, is 
therefore negative. The magnetically driven angular momen- 
tum transport is thus outward in both cases. 

Unlike the AMRI, the TI yields right-handed spirals (Fig. 
fTTT i. The pattern in the bottom row (J3 — 10) has an az- 
imuthal wavenumber m = 1, in accordance with the insta- 
bility map Fig. [8] The top row, however, represents a pattern 
with m = 2, which also exists in the nonlinear regime as 
predicted by the bottom plot of Fig. [8] In this case there is 
no clear correlation between the radial and azimuthal compo- 
nents of the field perturbations. 

The nonlinear simulations, of course, do also provide the 
amplitudes of the fields in the resulting magnetic pattern. Here 
we only note the overall result that the AMRI produces much 
higher field strengths than the TI. One might speculate that the 
AMRI exists due to the differential rotation which is always 
able to induce strong fields but a detailed study of the energy 
aspects of the magnetic instabilities is out of the scope of the 
present paper. 
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FIG. 1 1 : Components of the magnetic pattern (left: radial component, right: azimuthal component) for TI (slow rotation). Re = 30, Ha = 130. 
Top: j3 = 2, Bottom: f3 = 10. The fields are normalized with Bi n . [Ib = 1, Pm = 1. 



